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Abstract. 

In recent years, PHENIX has studied many important observables related to heavy- 
flavor physics through their leptonic decay measurements including the invariant yield 
of electrons from nonphotonic sources, and prompt single muons, both of which are 
dominated by D and B mesons. Charm and beauty cross-sections were measured and 
compared through single lepton, and lepton-hadron correlations in p+p collisions at 
•\/s = 200 GeV. Observables for quarkonia production such as invariant yield and 
polarization were also measured in p+p collisions. In Au-|-Au collisions, preliminary 
results for the Raa for single electrons and a 90% CL upper limit for the suppression 
of Ts were produced. And in d-fAu collisions, a preliminary Rcp study for J/ip 
production in different centrality ranges was extracted. 
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1. Introduction 

Measurements of heavy-flavor production in p+p interactions at collider energies serve 
as important tests for perturbative quantum chromodynamics (pQCD). It also provides 
an important baseline for study of the medium created in relativistic heavy ion collisions, 
since the transition from hadronic matter to a quark gluon plasma at high temperature 
is expected to be achieved in such collisions. On the other hand, the hadronization 
of partons in vacuum or cold nuclear matter, represents a non perturbative QCD 
process. In order to measure the baseline, disentangle the cold nuclear matter effects and 
determine the hot nuclear matter effects, the Relativistic Heavy Ion Collider (RHIC) 
provides, respectively, the colliding species: p+p, d+Kn and Au-|-Au or Cu-|-Cu. 

The PHENIX [1] experiment at RHIC has capabilities to measure both the leptonic 
and hadronic decay modes of the particles produced in the collisions. That allows 
a comprehensive study of open heavy-flavor and heavy quarkonia. Electrons are 
identifled in two central arms spectrometers covering the pseudorapidity range \ri\ < 
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0.35 and 2 X7r/2 in azimuth. Drift chambers and pad chambers are used for charged 
particle tracking, a ring imaging Cerenkov detector for electron identification and an 
electromagnetic calorimeter for electron identification and trigger. Muons are measured 
in two separate spectrometers with full azimuthal coverage and in the pseudorapidity 
range 1.2< \rj\ < 2.4. They are identified in panels formed by larocci tubes placed 
in front of steel plates for hadron absorption. While the muon arm absorber removes 
most of the hadrons, they are also measured in the central arms. In addition, beam- 
beam counters positioned at pseudorapidity 3.1 < \ri\ < 3.9, measure the position of the 
collision vertex along the beam direction and provide the interaction trigger. In A+A 
collisions, they also provide the event centrality characterization through a correlation 
of the beam-beam counters charge sum and the total energy deposited in the two zero 
degree hadronic calorimeters, located at ~ 18 m from the interaction point. 

2. Open heavy-flavor 

In PHENIX, since D or B mesons are not fully reconstructed, their yield can be 
statiscally derived by their semileptonic decays contributing to lepton spectra. These 
indirect measurements serve as baseline for both charm and bottom production in d+Au 
and A+A collisions. 

2.1. Single leptons 

The spectrum of electrons from heavy-flavor decay in p+p collisions was revisited [2] 
as an additional contribution to the previous spectrum of electrons [3]. The inclusive 
electron spectra consist of nonphotonic electrons from heavy-flavor decays; photonic 
background from Dalitz decays and photon conversions in the beam pipe; nonphotonic 
background from K^s and dielectron decays of vector mesons. The Kg^ contribution is 
small compared to the photonic background from Dalitz decays. The electron spectrum 
from heavy-flavor decay was determined by subtracting the background components 
from the inclusive spectrum using two independent and complementary methods. The 
electron spectrum from all known sources except semileptonic decay of heavy-flavor is 
calculated using a Monte Carlo simulation and subtracted from the inclusive spectrum 
in the cocktail method. The dominant source of background is the vr" Dalitz decay, so 
the measured 7r° and tt"^ spectra are used as input to the generator. The cocktail also 
includes contribution from quarkonium {J/ip, T) and the Drell-Yan process, which were 
neglected in |3lll] . Figure [1] shows the nonphotonic electron yield after the subtraction 
of the J/ip (that accounts for up to 16% of the electron signal at > 5 GeV/c), the T 
and the Drell-Yan contributions (both found to be negligible) in p+p collisions [2] . The 
bottom part of the figure shows the ratio of data to FONLL calculations [5]. From the 
nonphotonic electron spectrum derived in [3j the total charm production cross-section 
was deduced to be dec = 567 ±57 (stat) ± 193 (sys) fib. 
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Open Heavy Flavor: Charmonlum Subtracted 




Similar measurements were performed 
for single muons in which a background 
cocktail approach is also used. Back- 
grounds in the muon arms are divided into 
muons that result from the weak decay 
of light hadrons before the first absorber 
material and hadrons that penetrate the 
steel absorber to reach the deepest layer 
in the muon arm. Muons from heavy- 
flavor meson decay originate < 1mm from 
the collision vertex, while yields of muons 
from hadron decay exhibit a linear colli- 
sion vertex dependence before the first ab- 
sorber material due to decay kinematics. 
The background is then subtracted using 
vertex distribution and hadrons which are 
stopped in the muon identifier. 

Heavy-flavor single muon measure- 
ments in p+p collisions at y/s = 200 GeV 
and for ( y ) = 1.65 were reported in [6], 
with a corresponding diferential cross 
section for charm quark production at 
forward rapidity of dacc/dy\y=,Q = 0.243 ± 
0.013 (stat) ± 0.105 (data syst.) ^^^ff 
(PYTHIA syst.) mb. The new result 
( y ) = 1.65 was obtained over the 
extended range of 1 < px < 5.0 GeV/c 
[7]. Single muons are measured in both 
forward and backward directions, with the 
two spectra being combined to form the 
spectrum shown in Figure [2] [7J, which is 
compared to a FONLL calculation over the 
same rapidity. 

In Au+Au collisions, the heavy-fiavor 
signal and the ratio of nonphotonic to pho- 
tonic electrons were determined via exactly 
the same independent and complementary 
methods used in [3]. Similarly to the p+p 
case, the contribution of the J/ip to the 
electron spectrum in Au+Au collisions was 

determined |2]. Two conservative estimations of the upper and the lower bound of the 
J/ ip yield at high px were made using the measured nuclear modification factor Raa |8] 
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Figure 1. (a) Invariant cross sec- 
tions of electrons after subtrac- 
tion of the J ^ ^ Drell Yan con- 
tributions [2^. (b) Ratio of the 
data and the FONLL calculation. 
The upper (lower) curve shows 
the theoretical upper (lower) limit 
of the FONLL calculation. 
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Figure 2. Invariant cross sec- 
tions of muons from heavy-flavor 
decay at (y) — 1.65 7 . (b) Ra- 
tio of the data and the FONLL 
calculation. The upper (lower) 
curve shows the theoretical upper 
(lower) limit of the FONLL calcu- 
lation. 
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Figure 3. Raa of nonphotonic 
electrons with and witliout J/ip 
contribution [5]. 
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in the 0-20 % centrality range. 
PHENIX has previously presented the Raa 
result in the 0-10 % centrality range |1] 
showing a suppression level for the heavy- 
flavor electrons almost the same as for 7r° 
and for r) in high pt region. For the 0-20 
% centrality range, as that how J/ip was 
measured, the result with and without its 
contribution can be seen in Figure 3. Both 
the p+p and the Au+Au single electron 
yields from heavy-flavor decays decreased 
after subtracting the J/ip contribution, 
resulting on a small change in the Raa- 

2.2. Diparticle correlations 

In PHENIX, D meson decays cannot be 
fully reconstructed because of the small 
signal/background ratio and the poor PID 
of kaons. The ratio of ih — )■ e) to [(c — )• 
e) + (6^e)] function of electron p-p in 
the range 2 < pt <7 GeV/c was extracted 
from the correlation between the heavy- 
flavor electrons and associated hadrons and 
was used to separate the beauty from the 
charm contribution [9]. The extraction 
was based on partial reconstruction of 
the D/D e^K^X decay. The 

electron-hadron invariant mass is formed 
and exhibits a clear peak below the D 
mass. The shape and the efficiency of this 
peak depend on the charm or the beauty 
origin of the electron-hadron pair. The 
beauty decays present a lower efficiency 
and broader shape, because of their mass 

and alowed range of hadrons. In Figure 4, the single electron spectra for charm and 
bottom were measured from the ratio (6 — )■ e)/[(c — )■ e) -|- (6 — )■ e)], and the spectrum 
from heavy-flavor decays. 

Electron-muon azimuthal correlations in p+p collisions have proved to be a new 
method for the cc measurement at PHENIX [lOj. The heavy-flavor signal consists of 
opposite sign electron-muon pairs with A0 ^ tt, product of back-to-back charm pairs 
produced mainly via gluon fusion. 
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Figure 4. Invariant cross sec- 
tions of electrons from charm and 
bottom [9] with the FONLL cal- 
culation [5]. 
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e'-ti*, e*-|i' Azimuthal Correlations 



There is no contribution from Drell- 
Yan, thermal production, resonance de- 
cays for the background. Like-sign subtrac- 
tion removes combinatorial background 
and electron-muon pairs from dijets. Sub- 
stantial backgrounds in analysis are due to 
misidentified muons: hadrons that punch 
through absorbers and muons from light 
meson decays rather than from heavy- 
flavor. The background from hadron de- 
cays into muons was estimated and re- 
moved from the A0 correlation distribu- 
tion. Figure 5 shows the background sub- 
tracted A0 distribution. The band is error from decay muon subtraction, boxes are 
errors from electron-hadron subtraction. This new method provided an invariant yield 
of 2.11 X 10^''± 3.4 X 10"^ (stat.) ± 3.5 x 10"^ (sys.) and is a proof of principle for 
future charm cross section measurement in intermediate rapidity. 
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Figure 5. Electron-muon corre- 
lations from heavy-flavor [lOj . 



3. Quarkonia 
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Figure 6. (a) Preliminary measurement of J/ip yield multiplied by the branching ratio 
versus rapidity with both COM [IT] and s-channel cut Color Singlet Model (CSM) [H] 
calculations, (b) pr dependence of J/^ and tp' at midrapidity, together with the J/tp 
to tp' cross sections ratio ^ISj compared to that from other experimental results [14l|15j 



Using the 2006 data taking period, PHENIX has recently released a set of 
measurements on J/ip production in p+p collisions [13J including new results on the 
Pt dependence of its decay angular distribution in different reference frames [16]. The 
record of an appoximate factor of three of integrated luminosity, compared to the 
published one [17], allowed not only to extract the differential cross section of J/ip, 
but also the first ip' Pt dependent result at RHIC energies [13]. Figure [6] (a) shows 
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the preliminary measurement of the J/ip invariant yield versus rapidity, where the 
uncertainties are dominated by systematics with both forward and midrapidity yields 
well described by the s-channel cut Color Singlet Model (CSM) [12]. Figure [6] (b) shows 
the differential cross section of both J/ip and ip'. The ip' to J/ip cross sections ratio is 
shown in the lower pannel and compared to that from other experimental results [T^[T5] . 
The feed-down contributions to J/ip from -0' is 8.6 ± 2.5% and from Xc is < 42% (90 % 
CL.) and both are in agreement with world average results |18j . 
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(a) (b) 
Figure 7. (a) J/ip polarization parameter versus transverse momentum in the HX 

frame and theoretical calculations. Boxes are correlated point-to-point systematic 

uncertainties, (b) J/ip polarization parameter versus transverse momentum in the GJ 

frame. Boxes are correlated point-to-point systematic uncertainties [16] . 



With the same improved data set, the measurement of the transverse momentum 
dependence of the inclusive J/ip polarization can help elucidate the production 
mechanisms [TS|. Moreover, it is expected that the polarization of J/ip is modified by 
the presence of nuclear matter effects in d+Au collisions and hot and dense matter in 
Au-|-Au collisions [12]. The polarization was studied in three reference frames: helicity 
(HX), Gottfried- Jackson (GJ) and CoUins-Soper (CS) for < 5 GeV/c and \y\ <0.35. 
Both results for the polarization parame- ^ 
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ter in the HX and GJ frames are consistent I 
with zero for the full pt range and with a '^t 
1.8 sigma trend towards longitudinal polar- < 
ization for the J/ip above 2 GeV/c. How- ' 
ever, in the CS frame, no conclusion can 
be drawn due to its limited acceptance, be- 
sides the uncertainties in the current data. 
Figure [7] (a) shows the transverse momen- 
tum dependence of the J/ip polarization 
parameter in the HX frame together with 

currently available theoretical models: COM [20] and the s-channel cut CSM [12] . 



0,! 0,* 0.6 0-8 1 1.2 1,4 1.6 1.B 2 

Figure 8. Probability distribu- 
tion oiRAA in the T(15+2S'+35) 
mass range. 
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Figure [7] (b) shows the transverse momen- 
tum dependence of the J/tp polarization 
parameter in the GJ frame. 

PHENIX has recently presented pre- 
liminary results of dielectrons in the 
T(l^ + 2S + 35) mass range of [8.5-11.5] 
GeV/c^ in p+p and Au+Au collisions. Fig- 
ure 8 shows the nuclear modification factor 
Raa which was evaluated with Poisson dis- 
tribution and was found to be lower than 
0.64 in 90% confidence level. This result 
may suggest that the T states 2S-I-3S are 
melted in the medium formed in Au-|-Au 
collisions. However, no final conclusion can 
be drawn before the knowledge of the cold 
nuclear matter effects in d+Au collisions 
and a better understanding of the contin- 
uum contribution in this mass range. 

During the 2008 data taking period 
with d+An collisions, PHENIX recorded 
a factor of 30 higher statistics than the 
2003 d+An Run. Preliminary studies of 
the J/ip yield ratio between central and 
peripheral events {Rcp), function of 

rapidity, show that the current parton modification functions and uniform breakup cross 
section in the hadronic matter cannot describe the behavior of the J/tp suppression in 
d+An collisions. Figure 9 shows the measured Rcp in three centrality ranges, together 
with the prediction assuming gluon distribuion modification EKS98 [21] for different 
breakup cross sections (a, 0-15 mb). 
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Figure 9. Central to periph- 
eral J/tjj production as a func- 
tion of rapidity in three centrality 
ranges, together with the predic- 
tion assuming gluon distribuion 
modification EKS98 [21] for dif- 
ferent breakup cross sections (cr, 
0-15 mb). 



4. Conclusions 



PHENIX has measured single electrons in both p+p and Au-|-Au collisions at 
midrapidity and single muons in p+p collisions at forward rapidity at a/stvtv = 200 
GeV. The ratio of the yield of electrons from bottom to that from charm in p+p 
collisions at the same value of the cm. energy provided the first measurement of the 
spectrum of electrons from bottom at RHIC. All these results were presented in the 
context of pQCD calculations. The electron-muon azimuthal correlation result provided 
a clear measurement of heavy flavor production in a rapidity range not yet studied. 
Quarkonia yields in p+p collisions included the J/ip, the T, and the first tp' differential 
cross sections at RHIC energies, besides the first J/ip polarization measurement for 
two polarization frames, providing additional insight into the understanding of the 
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production mechanisms. 
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